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ABSTRACT 

We systematically investigate the error sources for high-precision astrometry from 
adaptive optics based near-infrared imaging data. We focus on the application in the 
crowded stellar field in the Galactic Center. We show that at the level of < 100 /xas a 
number of effects are limiting the accuracy. Most important are the imperfectly sub- 
tracted seeing halos of neighboring stars, residual image distortions and unrecognized 
confusion of the target source with fainter sources in the background. Further contrib- 
utors to the error budget are the uncertainty in estimating the point spread function, 
the signal-to-noise ratio induced statistical uncertainty, coordinate transformation er- 
rors, the chromaticity of refraction in Earth's atmosphere, the post adaptive optics 
differential tilt jitter and anisoplanatism. For stars as bright as tok — 14, residual im- 
age distortions limit the astrometry, for fainter stars the limitation is set by the seeing 
halos of the surrounding stars. In order to improve the astrometry substantially at the 
current generation of telescopes, an adaptive optics system with high performance and 
weak seeing halos over a relatively small field (r < 3") is suited best. Furthermore, 
techniques to estimate or reconstruct the seeing halo could be promising. 

Key words: astrometry, infrared: stars, Galaxy: centre, instrumentation: adaptive 
optics, techniques: high-angular resolution 



1 INTRODUCTION 

The Galactic Center (GC) is a unique celestial laboratory. 
It hosts a massive black hole (MBH) the existence of which 
is now widely accepted in astronomy (WoUman et al.||1977 



[Genzel et aL|l996 , Ghez et al.|1998 : Schodel ct al.|2002| ). The 



tightest constraints on the mass have been put by observing 



the motions of individual stars as they orbit the MBH ( Ghez 
|et al.|2008||Gillessen et aL|2009[ ). The statistical error of the 
mass is as low as ~ 1.5%, the systematic error is considerably 
higher with ~ 10% owed to the uncertainty in 7?o of ~ 5%, 
the distance to the GC. 

These astonishing measurements rely on high-angular 
resolution, adaptive optics (AO) assisted, near-infrared ob- 
servations, obtained at large telescopes; namely the Keck 
and VLT facilities. From imaging data astrometric positions 
are derived, reaching an accuracy of ~ 300 /ias. This a fac- 
tor of 200 smaller than the maximum possible image resolu- 
tion of ~ 60 mas in K-band, while the apparent orbit of the 



most important star of these so-called S-stars, S2, measures 
~ 190 mas. 

With the black hole paradigm now being well estab- 
lished, the efforts turn towards detecting post-Newtonian 
effects. From General Relativity, one expects a prograde 
periastron shift of 0.2° per revolution of 15.8 years for S2 
(Ru bilar fc Eckart] |2001). This corresponds to an apparent 
position shift of « 800 /^as, which however is not as easily 
detected as it seems, since it needs to be measured from 
the same data from which also the orbital elements have to 
be determined. Secondly, a population of dark stellar rem- 
nants such as stellar mass black holes and neutron stars 
might be present around the MBH ( Morris|[l993 1. Such an 
extended mass component will lead to a retrograde perias- 
tron shift. If for example « 0.1% of the mass inside the 
S2 orbit were extended, the resulting retrograde shift would 
cancel the relativistic prograde precession ( Rubilar & Eckart 



2001 1. Detecting either of these shifts would be astrophysi- 



cally extremely interesting. 

Future ground-based telescopes will further increase the 
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spatial resolution accessible, thus potentially increasing the 
astrometric accuracy further. Trippe et al. (in prep.) investi- 
gate in general the astrometric performance of an extremely 
large telescope. Weinberg et al. ( 2005 1 simulated the ex- 



pected advance in monitoring stellar orbits in the GC when 
using an extremely large telescope, concluding that the up- 
coming facilities indeed will be capable of measuring the 
post-Newtonian effects and that mass of and distance to 
Sgr A* can be determined with unprecedented precision. 

High-precision astrometry is important for other sci- 
entific cases, too. Intermediate mass black holes might be 
found by detecting accelerations of stars in the vicinity of 
the black hole. Potential sites are for example lo Cen ( |Noyola| 



et al.|2008 Anderson fc van der Marel|2009J or the compact 



Schodel 



)09) 

stellar grou p GC IRS13E ( [Maillard et aIT ( |2004[ ) 

|et al.| ( [2"005[ ); Fritz et al. (in prep.)). Another area where AO- 
assisted high-resolution imaging profits from excellent astro- 
metric capabilities is the domain of binaries and substellar 



companions (|Chauvin et al.|(|2004||; 


Neuhauser et al.|(|2007|); 


Kohler et al. ( 2008 1; Lagrange et al. 


(20081; Neuhauser et al. 


(20081; Kervella et al. (in prep.)). 



This paper aims at investigating the limits of astrome- 
try, in particular in the GC. This is of great value both for 
current and future studies of stellar orbits. Our approach is 
to estimate the various error sources mostly empirically, i.e. 
from existing VLT data. The main parameters characteriz- 
ing the reachable accuracy for any given star are its magni- 
tude and distance from Sgr A*, and for any given image the 
Strehl ratio. The empirical numbers used were determined 
for the VLT NIR-AO imager NACO, operated in Kg-band 
with a pixel scale of 13mas/pix. The field of view being of 
interest here is Ofoi = ±2", given by the apparent size of the 
GC stellar system. We concentrate on the data set from 13 
March 2008, which was of good quality (Strehl ratio larger 
than 30%) and therefore is well-suited to study systematic 
effects^ In that sense, this paper is a case study, but the 
dependencies on the Strehl ratio which we exploit make the 
results applicable for a wider range of data. 



2 STATISTICAL UNCERTAINTY 

The fundamental limit to astrometry for diffraction limited 
data is given by the number of photons recorded and the 
image resolution, with higher SNR and smaller point spread 
function (PSF, FWHM O) increasing the positional accu- 
racy. For a circular aperture of size D and an observation 
wavelength A one has 

(1) 



Ofwhm ~ 1.028 



D 



Lindegren ( 1978 1 gives the following formula for the position 



= i - — . (2) 
TT D SNR ^ ' 

We have verified this formula explicitly by simulating dif- 
ferent SNR for a simple Gaussian PSF. The sampling was 
chosen to be similar to VLT data from imaging GC observa- 
tions with OpwHM ~ 6 pix which occurs for a Strehl ratio of 



We are using data from ESO program 179.B-0261. 
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Figure 1. Left: SNR for a source with FWHM = 6 pix and Gaus- 
sian shape in a noise image as a function of the integrating radius 
used, showing that a source area with r = 4 pix is optimal. Right: 
The data points show the positional uncertainty in simulated im- 
ages as function of the SNR. The solid line is the prediction from 
equation [2] 

18?([^ We generated a set of noise images, in each of which 
the PSF was added. First, we searched the optimum radius 
for the source region for the given sampling. We determined 
the SNR as ratio of fiux inside a certain radius and the stan- 
dard deviation of the fiux inside the same area in the noise 
image. The optimum was 4 pixels (figure [T] left), indepen- 
dent from the peak intensity. Next, we varied the SNR by 
adding the PSF several times to each noise map with the 
peak intensity ranging from 0.3 x to 200 x the noise level. 
We determined the position of the PSF by a Gaussian fit to 
each image. The standard deviation over the set of noise im- 
ages per intensity then measures the positional uncertainty 
for that intensity level. The corresponding SNR was calcu- 
lated with the optimum radius. We plot the positional error 
as a function of SNR in figure [ij right, which shows that our 
explicit simulation matches equation |2] While the 1/SNR 
behavior is obvious, this checks mainly that the numerical 
factor in equation [2] correctly describes the accuracy for the 
optimum signal extraction. 

Empirically, we find for the VLT data that a source in 
K-band creates a total current (photo electrons per second) 
on the detector (in its high dynamic mode with the readout 
scheme double read-reset-read) of 

-,0.4 (25.4-mK) 



SpixTie, obj [el/ sec] 



10" 



(3) 



If source photon noise were the only error source, for a 
TTiK = 14 source (i.e. S2), within one second the position 
uncertainty would drop to 50 /las. However, also the noise of 
the sky background and the read noise of the detector have 
to be taken into account in equation [2] 

In a single pixel in a single-pointing, combined frame 
the noise after the usual data reduction steps of sky sub- 
traction and fiat fielding should be (assuming that the flat 
field on average equals 1 and that A^sky frames with the same 
exposure time as the A^obj object frames were averaged for 
the sky) 
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^c, obj 



bj 
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obj 



sky 



X (fRON + J^c.sky), 

(4) 



^ This is taken from the actual data, spanning the years 2002- 
2009. Since 2007 we are getting consistently smaller FWHM val- 
ues and higher Strehl ratios. 
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Figure 2. The solid lines (referring to the left ordinate) show the 
noise contributions per pixel as function of illumination according 
to equation [4] for a single image, which actually is the combina- 
tion of two object and two sky frames of 17.2 sec each. (Our data 
is mostly obtained in a mode with 2 detector reads of 17.2 sec 
per data file, a setting that is chosen to avoid saturation effects 
for most of the stars and minimizes overheads.) The total noise 
at faint magnitudes is dominated by the read noise. For realistic 
data sets, we use more than two sky frames, such that the read 
noise contribution from the sky is reduced compared to this il- 
lustration. For brighter objects, the photon noise is the dominant 
contribution. The regime in which the noise of the flat becomes 
important already is in the highly non-linear regime of the detec- 
tor which normally is not used. For the conversion to magnitudes, 
a sampling of 6 pix per FWHM and a Gaussian PSF that contains 
30% of the light were assumed. The dashed line (referring to the 
right ordinate) gives the resulting positional error according to 
equation [2] for the total noise. 



where at — 0.15% is the empirically determined relative 
noise of the flat field, s = (rie.obj — ne.sky)// the reduced 
signal per pixel and cr ron = 46.2 the read noise of the de- 
tector in electrons. Typically, the sky brightness per pixel 
for the given setup is 



''c, sky 



[el/sec] =21.5 



(5) 



The solid lines in figure [2] illustrate the noise behavior 
for a single pixel as given by equation |4] With the noise 
(equation [4| and the signal (equation [3| , it is straight for- 
ward to evaluate equation [2] if one assumes a certain PSF 
shape. As a simple example, we show the resulting posi- 
tional uncertainty for a case in which the position is esti- 
mated from a Gaussian shaped PSF-core with a FWHM of 
6 pixels that contains 30% of the stellar light. As above, 
the hght from inside a 4-pixel radius was used. The error 
so obtained represents the statistical limit to the positional 
accuracy for a single, reduced frame. Usually, astrometry 
is done on combined objects frames with varying pointing 
positions ('mosaics'). This improves the statistical precision 
limit by a factor ^ Appointing compared to equation | 



3 ASTROMETRY IN THE GC 

We briefly describe the method by which we obtain astro- 
metric positions. The lack of any extragalactic background 
source in the NIR is one of the main obstacles for astrometry 
in the GC. Therefore all position measurements are relative 



to other sources in any given image. The link to the inter- 
national celestial reference frame ICRF is only possible due 
to a set of SiO maser stars, which are both NIR and radio 
sources. The position vectors relative to Sgr A* of the latter 



can be measured with high accuracy ( Reid et al. 2003 2007 1. 
The position and motion of Sgr A* in the ICRF in turn is 
well known ( |Reid et al.|[7999l [Reid fc Brunthaierl[2004 1 . 

An additional complication is that with the cur- 
rent Ikxlk NACO detector, a suitable sampling (~ 
4pix/FWHM= 13mas/pix, Trippe et al. in preparation) is 
only reached for the central few arcseconds (depending on 
the dithering scheme), while the SiO maser stars are found 
out to 20". Therefore it is more practical to measure the 
SiO maser positions with a larger pixel scale (27mas/pix for 
NACO) and relate the finer scale to the coarser sampling by 
a set of reference stars that can be reliably detected in both 
scales. That also helps to overcome the large dynamic range 
needed, given that the brightest SiO maser star used for as- 
trometry has a magnitude of ttik ~ 8.5 (Blum et al.||2003 



Reid et al.|2007 \ and that one is interested in the positions 
of stars as faint as niK ~ 19- 

In practical terms, the procedure is as follows: From 
a set of images, obtained between 2002 and 2009 in the 
27mas/pix scale, we derive astrometric positions and proper 
motions for a set of ~ 100 reference stars. This relies on the 



work of Reid et al. (20071, which allows us to calculate the 



astrometric positions of the SiO maser stars for the given 
NIR epochs. We use mosaics corrected for their geometric 



distortion ( Trippe et al. 2008 1 and a full, 6-parameter lin- 
ear transformation to relate pixel and astrometric positions 
of eight SiO masers. For any given image in the 13mas/pix 
scale we determine the PSF from the image and deconvolve 
it using the Lucy- Richardson algorithm ( Lucy|[T974 1 . After 
beam restoration with a Gaussian beam we determine the 
stellar pixel positions by Gaussian fits to the positions, both 
to the reference stars and the sources targeted. The transfor- 
mation used to link the astrometric reference star positions 
to their pixel positions is a 20-parameter, third order poly- 
nomial transformation, which should also implicitly correct 
for any large scale (5") image distortion. This is useful since 
for the smaller pixel scale we were not able to construct a 
reliable distortion model, indicative that the effect is fairly 
small (Trippe et al. 112008 1. For a more complete description 



of the procedure see Gillessen et al. (20091. 



4 SYSTEMATIC UNCERTAINTIES 

A multitude of systematic uncertainties are present in the as- 
trometric data, the most important ones being atmospheric 
turbulence, image distortions, unrecognized source confu- 
sion and uncertain PSF halos. The time scales involved are 
very different; atmospheric effects are present even in a sin- 
gle frame, while confusion of sources happens on time scales 
of years. The section discusses the error sources (roughly) 
by increasing time scale involved. 
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4.1 Fast atmospheric Limitations 

4.1.1 Anisoplanatism 

Time variable refraction in Eartli's atmosphere, induced by 
turbulence cells, blurs ground-based astronomical images, an 
effect whicii is called 'seeing'. For the observations discussed 
here, the seeing is partly corrected by the adaptive optics 
(AO) system. The resulting PSF is a superposition of a close 
to diffraction-limited core with a seeing-limited halo. Useful 
parameters to describe the performance of the AO are the 
Strehl ratio (SR; the ratio of measured central flux compared 
to the diffraction-limited central flux) and the FWHM of the 
PSF. The AO correction in our data is achieved with a single 
guide star (single-conjugate AO). 

The tip-tilt component of the wavefront errors between 
any two objects in the fleld of view gets less and less corre- 
lated with increasing distance between the two objects, an 
effect which is called anisoplanatism. Complete decorrela- 
tion is reached for the so-called isoplanatic angle, which has 
a typical value of 10" — 20". As a result of anisoplanatims, 
one observes random variations in the relative positions of 
any two objects. Anisoplanatism is seen as PSF elongation 
during individual exposures, and as residual, differential tilt 
jitter between successive frames. 




Distance {" 



Figure 3. Differential tilt jitter in the data from 13 March 2008. 
The dashed lines show linear fits to the positional scatter be- 
tween pairs of stars in the direction parallel or perpendicular to 
difference vector between the two stars. The empirical model used 
was = Uq + P^d^. The fitted slopes /3 for large distances are 
0.114 ± 0.002 mas/" (blue) and 0.071 ± 0.003 mas/" (red). The 
solid, black line and corresponding data show the remaining as- 
trometric effect in R.A. and Dec. after a full linear transformation 
(the range is smaller since positions are counted from the center 
of the field then). 



4.1.2 Differential Tilt Jitter 

The differential tilt jitter between any two stars in good 
approximation linearly increases with distance between the 
two stars. The effect averages out with exposure time. For a 
Komolgorov turbulence spectrum the dominant term of the 
differential tilt jitter scales with telescope diameter D and 



integration time t like (Cameron et al.|2009l 



„2 

'^±, TJ 



-7/3 

t 



(6) 



The effect is big ger in the direction connecting the two 
stars than perpendicular to it. The time constant r char- 
acterizes the time it takes a turbulence cell to move over 
the telescope aperture. The constant a is related to the sec- 
ond moment of the atmospheric turbulence profile. A typical 



value, following from the numbers in Cameron et al. (20091, 
is a « 3 when a is measured in mas, D in meters and 8 in 
arcseconds. 

We looked for the effect in the data set from 13 March 
2008 by calculating the variance of stellar distances in a 
sequence of subsequent exposures that were obtained at ex- 
actly the same pointing position. We decomposed the differ- 
ence vectors into the directions parallel and perpendicular 
to the vector connecting the two respective stars. Figure [3] 
(dashed lines) shows the resulting scatter as a function of 
distance between the two stars. Clearly the differential tilt 
jitter is visible and the measured ratio cry /(jj_ = 1.62 ± 0.05 
is consistent with the expected value of \/3- Since we do not 
have the data to evaluate the turbulence profile for our ob- 
servation, we can check for plausibility only. From the slope 
of the relation doTi/dO — 0.071 ± 0.03 and given the expo- 
sure time of 34.4 s we estimate the wind crossing time over 
the telescope aperture of 8 m to be 0.4 s, corresponding to a 
wind speed of u ~ 20m/s, which seems reasonable. 



Since the field of interest OfoI = ±2" is practically al- 
ways smaller than the isoplanatic angle, the tilt jitter will 
act in a correlated way on all of the stars. Hence, the ef- 
fect will to a large extent cancel out for astrometric posi- 
tions defined as relative positions to a set of reference stars 
(Cameron ct al. 2009). The effect of the corresponding full 
linear transformation (i.e. including shear terms) is shown 
for the test data in figure [s] (solid lines). The remaining as- 
trometric error for a single frame is < 0.2 mas and in the 
center even < 0.1 mas. Due to the mosaicking, that number 
will decrease further by \/ Nf^^^^s and for a typical data set 
of > 70 images differential tilt jitter should contribute less 
than 20/ias. Hence, using a simple linear transformation is 
sufficient to effectively eliminate this error source. The more 



elaborate approach by Cameron et al. ( 2009 1 is not neces- 



sary for our current GC data sets. 



4.1.3 PSF elongation 

The AO corrects the tip-tilt jitter optimally for the the guide 
star. Since the exposure times are much longer than the time 
constant for the AO, the image retrieved can be thought as 
a superposition of many short images, each of which is cen- 
tered on the guide star. As a result, the random errors from 
equation[6]vanish at the guide star and the image is sharpest 
there. Further away, the integrated jitter increases the PSF 
width with distance to the guide star, and the particular 
form of equation [6] leads to an elongation of the PSF in 
direction towards the guide star. Since the analysis uses a 
constant PSF (either for deconvolution or PSF fitting) when 
estimating stellar positions, the effect will lead to increased 
position errors. 

We assessed the error by comparing the differences of 
extracted position from various PSF estimates, obtained 
from stellar images with similar brightnesses at different po- 
sitions in the field of view. We used the high-quality data 
set from 13 March 2008. Using two different PSF estimates 
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Figure 4. Empirical position errors as a function of stellar magni- 
tude, estimated using a PSF extracted from the respective data. 
The red, dashed line indicates the errors for all stars found by 
starfinder in a single frame and is dominated by anisoplanatism 
at brighter magnitudes. The red, solid line shows the errors for 
isolated stars obtained from starfinder on the same single image in 
a restricted area (the central 2", Spoi). The blue line indicates for 
the same stars the errors obtained by applying a Lucy-Richardson 
deconvolution to the single frame, and fitting the sources with 
Gaussian functions. The black, dashed line shows the errors ob- 
tained by simulating isolated stars and finding them back. The 
black dotted line is showing the expectation for a mosaic, obtained 
by describing the blue line with an empirical function and square 
root scaling with the number of frames used in the mosaic. The 
black, solid line shows the measured errors on the mosaic which 
flatten at a level of 70 fias. The data points show the measured, 
statistical errors for the mosaic from 13 March 2008. 



from the same region leads to a typical position difference of 
50 /xas. When the two PSF estimates are extracted from two 
samples 3.5" apart, we only find a small increase of the typ- 
ical position difference to 70 ^as. At 6" distance, the effect 
starts to hurt more severely since it reaches ~ 180 /las. At 
7" distance, it even reaches 400 fias. These values illustrate 
the effect of the isoplanatic angle. It is worth noting that 
anisoplanatism can be significantly more severe for worse 
atmospheric conditions. 

For the analysis, the effect will contribute only at a 
very low level, if the set of stars from which the PSF is con- 
structed is chosen carefully; i.e. within QfoI ~ ±2" of the 
target region. In GC data sets, this can easily be achieved 
and therefore we expect anisoplanatism to be only a very 
small error contribution. If the size of the target region ex- 



ceeds 4" , as for example in the work of Trippe et al. ( 2008 1 
or |Schodel et al.] ( |2009[ ), the anisoplanatic effects should be 
taken into account for astrometric measurements by varying 
the PSF over the field. 



4-1-4 Comparison with data 

Figure |4] compares the expected position errors with the ac- 
tual scatter in the data, using again the data set from 13 
March 2008. We obtained the expected errors by simulating 
isolated stars and remeasuring their positions, using a noise 
floor as measured in single frames and using the respective 
PSF also obtained from the data. The position scatter be- 
tween subsequent single frames is a useful estimator for the 
position error due to the statistical limit and atmospheric 
turbulence, since the frames were taken with exactly the 



same pointing position. This means neither image distor- 
tions nor confusion events bias the scatter. 

We find that compared to the simulations, positions 
found by starfinder (Diolaiti et al.|2000[ ) scatter as expected 
from anisoplanatism. The scatter gets reduced, if one re- 
stricts the analysis to stars that are sufficiently isolated in 
a small area {Qfoi = ±2"). For ttik > 15, the scatter is 
only moderately higher than in the simulation. In compari- 
son, our technique of Lucy-deconvolving the image performs 
a bit better than starfinder and is for rriK > 14 essentially 
as good as predicted by the simulation. The reason why 
starfinder is performing a bit worse might be that it was 
allowed to calculate a (spatially varying) background map 
for each single image, which however is not needed for the 
data. The additional degree of freedom then introduces ex- 
tra noise, from which the deconvolution technique does not 
suffer. Both for starfinder and deconvolution, the brighter 
stars have a higher scatter than what one would expect 
from the simulations and the values reach a floor. Its level 
of « 125 ^as is consistent with being the combination of 



PSF uncertainty (section 4.4.1 1 and differential tilt jitter as 
present in figure |3] (data after transformation). 

The empirical curve for the deconvolution of a single 
frame was then parametrized by a simple model of type 

= aa^{mK) + (^o- For the given data, a = 1.1 and 
(Jo = 100 fias are an excellent description of the measured 
errors. We scaled this relation by \/ A''frames = \/24 in order 
to estimate the expected positional uncertainty for a single- 
pointing mosaic (figure |4|. For rriK > 15 the actual data 
(obtained again with deconvolution) follows the expected 
relation and flattens for sources with ttik < 14 at about 
70 fias. This is the value as expected from the PSF uncer- 



tainty (section 4.4.11 and the anisoplanatism over the fleld. 
Hence, we conclude that no additional terms come in by 
combining multiple frames to a single-pointing mosaic. 

For comparison, we also show in figure [4] (data points) 
the positional errors as obtained by comparing two mosaics 
created from subsamples using the data from 13 March 2008. 
The match with the simulated is actually a combination of 
the fact that 1.5x more single frames were used here for each 
submosaic, and that the error is composed of the fitting error 
plus the positional difference between the two submosaics. 



4.2 Slow atmospheric Limitations 

4-2-1 Differential refraction 

Refraction in Earth's atmosphere displaces stellar images 
towards the zenith by 



1 pTo 



2nl poT 



tan 2 ~ 44 tan z 



(7) 



where z is the zenith angle of the observation and the con- 
stant of 44" is the refraction constant for typical conditions 
at Paranal as evaluated from the refractive index no at stan- 
dard conditions of po — 1013.25 mbar and To — 273.15 K, 
scaled by the current pressure and temperature. Since we 
are using relative astrometry, only the differential effect is 
relevant. It amounts to 



SR ^ eFoi[rad] x 44"/cos^ z 



(8) 



Numerically, this evaluates to values around 3.5 mas for 
z = 60° and 0poi ~ ±2". It would thus be a large error 
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Figure 5. Assuming black body emission, the astrometric shift of 
a source due to the chromatic effects of Earth's atmosphere com- 
pared to a fictitious emitter of T = 10^ K is shown as a function 
of temperature. The red lines assume a K-band measurement, the 
blue lines H-band. The solid lines are for a zenith angle of 20°, 
the dashed lines for z = 40° . 
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Figure 6. Positional shift of a black body source per 0.1 mag 
change in extinction as a function of zenith angle as it results from 
the chromaticity of Earth's atmosphere. The red lines show that 
the effect for K-band data is nearly independent from the source 
temperature. For H-band data (blue lines), the effect slightly in- 
creases with temperature. 



source. However, the effect is absorbed into the shear terms 
of the transformation; a full linear transformation leaves a 
remaining second order effect that is less than 3 ^as. 

4-2.2 Chromatic effects 

The dependence of the refractive index n on the observation 



wavelength (e.g. as given by Edlen (1953l) for the H- and 
K-band and typical Parana! conditions is approximated by 

10^ = 2029.94-2.87A2 + 2.16A2 



ni ■ 



I.44A2 + O.92A2 (9) 



where ni = n — 1 and A2 = A[/^m] — 2. Thus the atmosphere 
acts as weak spectrometer, as a result of which the effective 
position for a source is the weighted average of the wave- 
length dependent positions over the band ( |Helminiak|2009[ ). 
The weighting factor is the number of photons per wave- 
length which in turn is the convolution of the input spectrum 
with the interstellar extinction and the atmospheric and in- 
strumental transmission numbers. The size of the net effect 
grows quadratically with the bandpass used and hence nar- 
row band filters help to suppress the atmospheric chromatic 
effects on the astrometry. 

Figure |5] shows the resulting effect for black body type 
emission compared to a fictitious object with 10^ K. For 
young, hot stars (most of the S-stars have T ~ 25000 K, 
Ghez et al.] ( |2003[ ); [Martins et"ar| ( |2008[ )) the K-band mea- 
sured positions for typical zenith angles are altered by 
~ 20 fias. For stars with a temperature of 6000 K the K-band 
shift would be ~ 70 /ias. Below 5000 K, the approximation 
as black body breaks down due to the presence of the broad 
CO absorption features around 2.3 fim, and the shift does 
not get larger anymore but rather is close to with a typi- 
cal spread of < ±20 fias compared to the fictitious emitter. 

Since our position measurements are relative to other 
sources in the field, the differential effect between reference 
stars and target sources is an error source for the positions 
obtained. The reference frame is constructed from early- and 
late-type stars in roughly equal numbers, such that in K- 
band it refers to coordinates shifted by « 10 pas compared to 
the fictitious object. Individual stars are then shifted again 
depending on their stellar type by ~ 10 /las compared to 



the coordinate system. Only G-type stars would experience 
a shift of ~ 50/xas; but main sequence stars of that type are 
by far too faint to be detected in the GC field. The net effect 
for K-band based astrometry is thus very small. 

For H-band data the effect is larger in the first place 
(figure |5| and secondly late-type stars actually can be ap- 
proximated by black bodies in H-band. For z — 40° and 
a late-type star the chromatic shift can reach 400 fias com- 
pared to the 10^ K emitter. For hot, young stars a more 
typical value is 50 fias. The astrometric net effect will thus 
be ~ 175 fias. 

Of course, these values strongly depend on the actual 
zenith angle and the exact input spectrum. Also, this error 
source is correctable, given that the zenith angle is known 
and if the input spectrum is at least approximately known. 

More subtle is the combined effect of the patchy na- 
ture of the extinction screen towards the GC and the chro- 
maticity of the atmosphere. Extinction variations lead to a 
change in observed color for each object as it moves behind 
the screen, the change in color in turn leads to a positional 
offset. By looking at the scatter of the observed K-band 
magnitudes of some isolated, bright stars (e.g. S8, SIO, S30, 
S65, S87 in the nomenclature of Gillessen et al. (20091) we 



conclude that the extinction variation AAk 0.1 for the 
central field of interest Opoi. Figure [6] plots the positional 
shift per 0.1 mag extinction variation as a function of zenith 
angle. For K-band data, the effect is ax < 20 pas, for H- 
band data ax < 100 pas. Also note that for larger regions in 
the GC of a few arcsecond, values of AE ~ 0.5 are reported 
( Buchholz et al.|2009 \ , which yields a correspondingly larger 
chromatic effect. 

Finally, it is worth investigating the chromatic effects 
for Sgr A*, the MBH itself, which in the NIR is a vari- 
able source powered by the synchrotron emission of rela- 
tivistic electrons and pos sibly changes its power law index 



P {vSu ~ v'^) with fiux (Eisenhauer et al 



Hornstein et al. 



|et al.|2006[|Krabbe et al. 2006[ ), but see also 
(20071. The astrometric effect ax per unit change of /3 is 
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Gillessen 



shown in figure [7] as a function of zenith angle. The effect 
reaches ~ 200 pas for a zenith angle of > 30° and a change 
of A/9 = 2. This is still a bit smaller than the currently 
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Figure 7. Positional shift in Kg-band of a power law type spec- 
tral energy distribution source (like Sgr A*) per change in power 
law slope A/3, uSi, ^ , as a function of zenith angle as induced 
by the chromatic nature of refraction in Earth's atmosphere. 
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Figure 8. SNR induced position error in a Kg-band mosaic of 96 
single frames as a function of Strehl ratio for stars of magnitudes 
between 11 and 18. The darker a star is, the more its SNR induced 
position error is sensitive to the atmospheric conditions. 



achieved astrometric accuracy, but it can be a big system- 
atic effect for future studies using extremely large telescopes. 
The easiest way to overcome the problem probably is the 
use of a narrow band filter, but also atmospheric dispersion 
compensators are feasible. Interferometric studies of Sgr A* 
( [Eisenhauer et al.,j2008j will use several spectral channels 
over the K-band. In each channel, the effective wavelength 
will be altered by a change in color of Sgr A*. Using five 
spectral pixels of A A = 0.1 /xm will already lower the as- 
trometric effect to < 10 /ias, and in addition a first order 
correction will be possible because the color of the source 
can be measured simultaneously. 



4-2.3 Influence of Strehl ratio 

For any given AO image, it is common practice to charac- 
terize its quality by the SR. It is strongly correlated with 
the FWHM of the PSF, and hence from equation ^ fol- 
lows that the SR limits the reachable accuracy. We sim- 
ulated the influence of the SR by placing measured PSFs 
into noisy background maps and finding these isolated stars 
back using starfinder. The resulting errors as a function of 
SR followed very well power laws with power law indices 
between —0.7 and —1.3 for stars between ttik = ff and 
niK = 18 (figure [sj . The power laws get steeper with in- 
creasing magnitude since the readout noise is more impor- 
tant for faint sources. As a consequence, the SNR induced 
position error for fainter stars is more sensitive to the ob- 
serving conditions than for brighter stars. As suggested by 
figure |8] and in practice, the SR induced position errors are 
smaller than other error terms, most notably the halo noise 
(section [4.4.4 1. Hence, our current data sets are not SNR 
limited. 



4.3 Instrumental Limitations 

4-3.1 Distortion 

Each imaging system suffers to some degree from optical 
distortions. These can either be a global effect for the image 
(such as classical optical aberrations) or small scale imper- 
fections. The global terms should mainly be absorbed by the 



third order transformation between pixel positions and as- 
trometric positions, but a floor of residuals is unavoidable. 
If we would use only one pointing position, these residuals 
would not influence the precision of the data, but would of 
course affect the accuracy. However, we use various pointing 
settings from run to run and even during one run (to increase 
area coverage). The residual image distortions effectively act 
like a random error. 

Our 13mas/pix observations are typically dithered in 
a scheme that four positions, offset by 3.5" from Sgr A* in 
both axes, are obtained. Since the field of view is ~ 13.5", 
the central « 7" are covered by each frame. From these 
data, we usually construct multi-pointing mosaics. The ef- 
fect of the image distortions here is a) that each source gets 
broadened slightly and b) that a superposition of the global 
distortions is present in the area which is covered by sev- 
eral pointings. Both the area of interest and the reference 
stars are taken from the central 7" which are covered by all 
pointings, and hence the resulting distortion pattern might 
be a complicated function, which we represent by a cubic 
transformation with 20 parameters. 

Following the approach in |Gillessen et ah] ( |2009[ ) we es- 
timate the size of the effect for the data from 13 March 2008 
by creating four single-pointing mosaics and comparing the 
distances of stars common in all pointings after having trans- 
formed them with a cubic transformation onto each other 
(figure [9|. The width of distribution of position differences 
then measures the residual image distortions, but includes 
also the effects of PSF uncertainty and differential tilt jitter. 
We obtained a value of = 148 /ias, of which « 50 /ias are 
due to the uncertainty in the PSF, ~ 40 /ias due to anisopla- 
natism, and ~ 125/\/24/ias « 25 /ias due to differential tilt 
jitter. The SNR induced position error is < 5 pias here, since 
the mean magnitude for the set of stars used was ruK = 13.1. 
Hence, we conclude that the residual image distortions affect 
the astrometry at a level of ~ 130 /ias. This value is smaller 
than the corresponding number in [Gillessen et al.| ( |2009[ ) 
who had not taken into account the other effects subtracted 
here and who were looking at a different data set from a 
time before the instrument got realigned. 

A different route to estimate the residual image distor- 
tions is to compare the high order terms of the transfor- 
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Figure 9. Compensation of image distortions by means of a cu- 
bic transformation for the data from 13 March 2008. The broad 
histogram is the distribution of differences of bona fide stellar 
positions from four pointing positions with a dither offset of 7" 
after having mapped the positions onto each other with a full lin- 
ear transformation. If the stars are mapped onto each other by a 
cubic transformation the same differences yield a much narrower 
distribution since image distortions are removed. The correspond- 
ing one-dimensional position errors ax are given, too. 



mations between several epochs. For this purpose we cal- 
culated the inverse of the third order transformations for 
each of the 73 NACO epochs with 13mas/pix data to con- 
vert a regular grid on sky back to the respective pixel po- 
sitions in each frame. Then we mapped these pixel posi- 
tions to a common epoch with a full linear transformation 
which was determined from the reference stars for the given 
epoch. Like this we take into account image shifts, rota- 
tions, scale changes and shear terms and therefore the re- 
sulting jitter of the grid points measures the variations of 
the higher order distortions. Since the set of reference stars 
on average still might be moving in a linear fashion we ac- 
tually determined the jitter as the residuals to a temporal 
linear fit for each grid point. We broke up the entire data 
set into three periods, corresponding to three hardware con- 
figurations: From 2002 to June 2004 a different detector was 
used than in the following years, and NACO underwent a 
major realignment late 2007. From the residuals we con- 
structed maps of the residual distortions (figure [l0| . For the 
first two periods, a typical value for the central arcsecond is 
CTj, « 300 ± 50 /ias, while for the last period a significant im- 
provement to ax ~ 120 ± 30 ^as is observed. This estimate 
agrees with the number quoted in the previous paragraph 
for the data set from 13 March 2008. 

For the 27mas/pix camera of NACO, the residual dis- 
tortions in absolute terms are higher. The procedure of 
[Trippe et al.| (2008) leads to residuals of order 1.2 mas, (see 
their figure 2, bottom; or figure 6 in Gillessen et al. (20091). 



Over time, the distortion appears stable, given that the dis- 
tortion parameters of [Trippe et al.| ( ["2008[ ) are consistent with 
each other for all data analyzed since 2002. For a distance 
of 10" from the optical axis (which differs from the center 
of the detector), the uncertainty of the geometric model is 
0.4 mas and hence is smaller than the residuals. We have not 
analyzed further, what limits the positional accuracy for the 
27mas/pix data. For our data sets, the residuals of this plate 
scale matter when the positions need to be related to an 
absolute coordinate system, for example for answering the 
question where radio-Sgr A* is located on the frames. For 
many questions, however, this does not matter, e.g. when 
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Figure 10. Maps of estimated residual image distortions. Con- 
tour labels are in mas. Top: For the period between June 2004 
and fall 2007. Bottom: For the period after the realignment in 
2007 until 2009. Note that identical color maps have been used, 
pronouncing the fact that the residual distortions are much less 
severe after the realignment. 



one is trying to detect accelerations. For a more detailed 



discussion see Gillessen et al. ( 2009 1 



4-3.2 Detector non-linearity 

Nominally, the current NACO detector is regarded as a lin- 
ear device up to 2/3 of its capacity, corresponding to 14000 
ADU for the readout mode of our data sets. Using a set of 
flat exposures with varying exposure times, we checked the 
linearity behavior of the detector and found that at 10000 
ADU the non-linearity is 7%, at 7000 ADU 4% and at 5000 
ADU still 2%. This would in turn mean that a PSF estimate 
is worse for brighter stars, possibly leading to positional bi- 
ases. Hence, we calculated a pixel-wise map of polynomial 
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Figure 11. Position error induced by the detector non-linearity 
under good observing conditions as a function of flux in ADU (top 
axis) and stellar magnitude (bottom axis) assuming an image with 
DIT= 17.2 s. The error bars show the standard deviation of the 
position errors and hence measure the scatter per bin. 



coefficients to correct the effect, which is now applied to all 
data. That correction doubles the usable dynamic range of 
the detector and thus allows the use of brighter stars. It is 
now linear to 2% up to 12 000 ADU, which for the data from 
13 March 2008 corresponds to ttik = 10.7. 

We also tested for the potential astrometric bias by 
comparing the positions from a linearity-corrected and an 
uncorrected frame. Figure [Tl] shows that the effect in a mo- 
saic for the brightest S-stars with ttik « 14 is < 50 /xas and 
decreases for fainter stars. 

Stars that are affected by saturation are offset sporadi- 
cally up to 100 X the usual positional error, i.e. in the multi- 
mas range. Such a big effect is easily detected in the data and 
hence such cases can be flagged as outliers. We thus don't 
assign a formal error here but rather expect that for some 
data sets, for bright stars saturation effects make the posi- 
tions unusable. Furthermore, heavily saturated stars (show- 
ing for example a flux plateau at the highest flux value) are 
not described at all by the PSF estimate; such stars are even 
recognized earlier, since the deconvolution does not yield a 
point-like structure in these cases. 



4.4 Analysis Limitations 

4-4-1 Knowlegde of the PSF 

Astrometry in a crowded stellar field necessarily needs as 
input an estimate for the PSF. Just using local maxima in 
the image as position estimates is not sufficient, since essen- 
tially all stars are located on the seeing halo of neighboring 
sources. This yields a background flux with a gradient and 
thus biases the positions of local maxima. Even worse, the 
background will vary with the AO correction from image 
to image and lead to a position jitter. Moreover, very faint 
stars in a seeing halo of a bright star do not produce a local 
maximum (but only change the shape of the seeing halo) and 
would be missed by an algorithm looking for local maxima. 

These problems are overcome by using the knowledge of 
the PSF, which can be obtained from the image itself with 



PSF fitting. A program like starflnder (Diolaiti et al 



20001 iteratively fits local maxima to the image until the 



image is sufficiently well represented by the set of sources 
found, each of which is characterized by a position and flux. 

• Deconvolution. The image on the detector is the con- 
volution of the PSF with the original object distribution. 
Since a convolution in image space corresponds to a multipli- 
cation in Fourier space, the inverse operation deconvolution 
in the Fourier domain is a division. The presence of noise 
and the limited range of spatial frequencies complicate the 
actual algorithm. Various methods are available and have 
been compared ( Ott et aLlllQQ Q). We mostly use the Lucy- 
Richardson algorithm ( |Lucyp974| ) for astrometry. 

Since both methods use the same input, one does not expect 
a major difference in astrometric performance. The fact that 
in figure |4] the starflnder errors are larger than the corre- 
sponding deconvolution errors might reflect a non-optimum 
user choice for the parameters of the starflnder algorithm 
and/or it might be due to the differences how the non-stellar 
background is dealt with. 

The fact that the input PSF is not perfectly known but 
has to be estimated needs to be accounted for in the error 
budget. In order to quantify the error, we used various sets 
of PSF stars and compared the resulting positions. Also, we 
compared the positions depending on the number of itera- 
tions used for the extraction of the PSF. From the scatter of 
the positions (after having transformed them linearly onto 
each other) we conclude that the uncertainty of the PSF 
knowlegde yields a positional error of SO/ias. 



4-i-2 Transformations 

The fact that we need to transform the measured pixel posi- 
tions with a set of reference stars is an error source because 
the transformation parameters also have errors. Fortunately, 
the corresponding statistical error can be lowered to a very 
low level. If A'^ stars are used as reference stars, each with 
a typical position error of ax, the additional error for the 
target sources due to the transformation is of order ax/\/N ■ 
By choosing a large enough sample of reference stars, this 
statistical error can thus be lowered to any desired value as 
long as a sufficient number of reference stars is available. Of 
course, this holds only if the reference stars can all be cho- 
sen from within the field in which the PSF can be treated as 
constant. Recently Gillessen et al. (20091 used A'^ « 100 ref- 



erence stars in the GC field, which suppresses effectively the 



an iterative procedure, e.g. using starfinder (Diolaiti et al 



transformation error to below 50 /^as. Fomer studies ( Schodel 
iet al._2002 ) had used A ~ 10 reference stars, which hence at 
the time was contributing at the few 10% level to the errors. 



4-4.3 Deconvolution 

We also investigated systematically the number of iterations 
to be used with our implementation of the Lucy-Richardson 
algorithm. Using a set of 66 well isolated stars in the data set 
from 20 July 2007 we obtained pixel positions for different 
deconvolution depths. 

For less than 500 iteration steps, the resulting positions 
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are biased up to 0.5 mas for bright stars compared to the 
positions in the undeconvolved image and compared to more 
deeply deconvolved images. For different stars the bias acts 
in different directions, and hence globally it appears to be a 
scatter of up to 0.5 mas. For 2000 steps, both bias and scatter 
for bright stars drop below 120 /las. The optimum is reached 
for ~ 10000 iterations; also fainter sources are unbiased with 
that number of iterations. Much larger values (like 50000 
steps) tend to put the flux of any star into a single pixel, 
which of course corresponds to larger positional errors again. 
The optimum depth in practice is between 5000 and 15000 
steps, with better data needing less deep deconvolution. 



4-4-4 Halo noise 

Visual inspection of the deconvolved frames shows that the 
background in the images is not flat but rather has a noise- 
like structure, see for example figure 1 in [Gillessen et al.| 
(20091. The deconvolution algorithm does not discriminate 



between sources and background and therefore any extended 
light present gets concentrated into this floor of faint local 
maxima. That light is due to the imperfections in the knowl- 
edge of the PSF wings (the dominant source) and real, phys- 
ical extended background light as expected from interstellar 
gas. Astrometrically, the spurious deconvolution peaks act 
as a noise source. It turns out that this halo noise is the 
main limitation for fainter sources. Experimentally indistin- 
guishable from this is the confusion noise due to unresolved 
background sources (section |4. 5. 1| |. 

While the term halo noise describes the influence of the 
imperfectly known PSF halo of surrounding stars, the PSF 
uncertainty discussed in section [4.4.1| refers to the effects of 
the imperfectly known PSF on the source itself. The PSF 
uncertainty from section [4 .4 . 1 1 would be present also fur suf- 
ficiently isolated sources, the importance of the halo noise is 
a consequence of the stellar crowding. It is also worth noting 
that the halo noise is not specific to the deconvolution tech- 
nique. Also a starfinder-like algorithm is not free from that 
noise source. It generally arises from the fact that the seeing 
halo of neighboring, brighter stars cannot be subtracted per- 
fectly, which owes to the large dynamic range in the images 
and the high surface density of stars. 

As a first test to estimate the effect of halo noise, we 
added Gaussian peaks with the same width as the typi- 
cal width for real stars into the deconvolved frame from 13 
March 2008 at random positions, avoiding positions at which 
the added source would overlap with an already existing and 
detected star. We then re-identified the random sources and 
determined their positional uncertainty by calculating the 
median deviation of the distribution of differences between 
the respective initial and re-identified positions. This tests 
by how much simulated stellar sources are shifted due to the 
combined deconvolution and confusion noise. The resulting 
positional error can be described very well by a relation of 
type 
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Figure 12. Position errors due to halo (and confusion) noise as a 
function of stellar magnitude and distance from Sgr A*. The er- 
rors are obtained by adding simulated stars to an undeconvolved 
frame and finding the sources back after deconvolution. The con- 
tour label units are mas. 



data set T.x has the following radial dependence: 
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where T,x{r) is a characteristic error for a source with niK ~ 
14 at a given radial distance r from Sgr A*. In the chosen 



The last bin, labelled M, corresponds to manually selected, 
well-behaved background regions in the deconvolved image. 
There are two reasons why the deconvolution and confusion 
noise decreases with radius: a) The unrecognized, faint stel- 
lar population falls off with radius, b) The stray light, lead- 
ing to a variable background on top of the PSF, follows the 
general trend that the light is roughly concentrated towards 
Sgr A*. 

In a second test we added the full PSF into the undecon- 
volved frame, neglecting the effects of anisoplanatism. The 
modified frame was then deconvolved as before and stars 
were identified, yielding again an estimate for the error in- 
troduced by the deconvolution with imperfectly known back- 
ground and PSF halo. This procedure yielded the errors as 
a function of magnitude and radius from Sgr A* as shown 
in figure |12[ The numbers obtained are broadly consistent 
with the results from the first test where Gaussian peaks 
were added to the deconvolved frames. 

Hence, the combination of halo and confusion noise is a 
severe limitation. Only for the brightest S-stars (mK ~ 14) 
it gets as small as the errors due to residual image distor- 
tions. In the next section, we show that the confusion noise 
actually contributes much less than the halo noise and we 
have to conclude that halo noise is limiting the astrometric 
accuracy for stars fainter than ttik ~ 15. 

Due to the importance of the halo noise, we have re- 
peated the last test (putting full PSFs into an existing im- 
age and finding the stars back) using PSF fitting (starfinder) 
instead of Lucy-deconvolution. We used 100 artificial stars 
per magnitude bin on the image from 13 March 2008, the 
positional error again is the width of the distribution of the 
differences between detected and simulated position. The 
results are very similar to what we found with deconvolu- 
tion. The positional uncertainty using starfinder is slightly 
larger than for the deconvolution, similar to the findings in 
figure |4] The photometric correctness on the other hand is 
slightly better for the PSF fitting method. This can be un- 
derstood, since the Lucy algorithm tends to pull fiux fro the 
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surroundings into local maxima during the iterative proce- 
dure. For a low SNR source this can bias the flux estimate. 
Overall, we conclude that halo noise is an important noise 
source, and it is not specific to the method used to extract 
positions. 



4.5 Astrophysical Limitations 

4.5.1 Confusion 

Astrometry in the GC is affected by the stellar crowding 
close to the MBH. The S-stars reside in a cusp with a steep 
density profile. As a consequence, the closer a star is to the 
MBH, the more frequently it is confused with some other 
source. Not all confusion events are recognized, and there- 
fore unrecognized confusion events will contribute to the po- 
sitional error budget. This error is already included in the 
discussion of section 14.4.41 It is nevertheless instructive to 
single out the effect of the stars, since that limit would re- 
main even if the instrumental halo noise could be removed. 

As an extreme example for a possible confusion error, it 
is worth mentioning that the orbit of the star S2 possibly was 
affected in 2002, during its pericenter passage, by such an 
event. Both recent analyses ( Ghez et al.|2008 Gillessen et al 



20091 therefore treat the respective 2002 data separately; 



either by ignoring it or by assigning large errors to it. 

In order to assess the magnitude of the confusion noise, 
we simulated stellar background populations in a Monte- 
Carlo fashion. That needed two basic input distributions: 
a K-band luminosity function and the radial surface den- 
sity profile. We based these on the findings of |Genzel et al.| 
(20031. We used three radial bins for our simulations: 0" < 
n < 0.2", 0.2" < ra < 0.8" and r-j. = 3.5". Since we had 
to assume the density also for stars much fainter than what 
can be measured, some assumptions had to be used. For 
ra we extrapolated the cluster KLF of |Genzel et al.| ( |2003[ ) 
down to rriK = 24. For ri,2 and ruK ^ 18 we used the KLF 
as estimated from the S-stars cusp, scaled to the respective 
expected surface density from the radial profile. At fainter 
magnitudes we extrapolated with a KLF from the same ra- 
dial region that only counts stars which are not identified as 
late-type stars. This essentially assumes that for rriK > 18 
only main sequence stars are present. Figure [13] shows the 
densities used. 

Using the assumed densities per magnitude we simu- 
lated stellar fields, using Gaussian profiles with a FWHM 
of 30 mas — 42 mas (magnitude dependent, to mimic the 
FWHM of real sources in the deconvolved frames) and a 
binning of 13mas/pix. For each image, the target star was 
placed in the center of a box of 12 pixels width and back- 
ground stars fainter than the target star were added. The 
positional error of the target star is given by the difference 
between input position and the position at which it is found 
back by fitting the star with a Gaussian profile plus a fioor. 
We used up to 10^ simulated images per magnitude bin. The 
distribution of positional differences per magnitude bin was 
then fit with a Gaussian, the width of which estimates the 
position error. The resulting errors as a function of magni- 
tude are shown in figure [Ts] 

For the S-stars cluster (r < 0.8") the error due to un- 
recognized confusion is of similar magnitude as the error due 
to halo noise. For larger radii, the confusion induced error 
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Figure 13. Assumed stellar densities and resulting position er- 
rors due to confusion noise: Dashed lines and right axis: assumed 
stellar surface density as a function of magnitude. Solid lines and 
left axis: resulting positional errors due to unrecognized confusion 
events with fainter stars. The blue lines are for 0" < r < 0.2", 
the red for 0.2" < r < 0.8" and the black for r = 3.5". 



is smaller than the halo noise as the stellar densities drop 
rapidly with radius. The increase of confusion error with 
stellar magnitude is well described by a power law of type 
C X 10°-^ ™K (as the halo noise). 



4-5.2 Relativistic Effects 

Currently, relativistic effects have not yet been detected in 
the data of any star orbiting the GC MBH. Actually, de- 
tecting effects will be possible probably first in radial 
velocity measurements ( Zucker et al.|2006 l. Astrometrically 
detectable deviations from Newton's law have not yet been 
seen in the GC. Vice versa, the effects can currently be ne- 
glected in the analysis. Still, a few effects are worth dis- 
cussing here. 



• Gravitational light deflection by the MBH. The 

deflection angle for a star sufficiently far behind the MBH 
(at distance z) can be approximated by 



6^20 ^ias - 





(11) 



where b is the impact parameter (Nusser & Broadhurst 



|2004[). F or smaller values oi z, 9 decreases compared to equa- 
tion |11| and for z = it is half of it. Generally the effect 
leads to very faint secondary images (rriK ~ 22 at best), and 
lensing only can become important when a star crosses the 
line of sight to Sgr A* behind Sgr A*. Currently none of 
the stars tracked is in that regime and we neglect the ef- 
fect. On the other hand, detecting the effects of secondary 
images would be extremely interesting. |Bozza fc Majicim] 



( 2005 2009 1 calculate the appearance of the secondary im- 
ages for the known orbits. It will be extremely difficult to 
detect these faint sources in the crowded field around Sgr A* 
ever. Another application would be the idea of |Alexander| 
(2001) who suggested to use the secondary images of back- 
ground stars to pinpoint the position of the MBH. 

• Gravitational light deflection by the Sun. Light 
deflection by the gravitational fleld of the sun can be ap- 
proximated by ( jLindegren" fc Bastian|2"006l ) 



61 « 4 mas x cot * /2 , 



(12) 
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where is the angle between Sun and the observed ob- 
ject. While the absolute magnitude of the effect is as large 
as 22 mas for 'I' = 20° (the minimum useful distance), the 
differential effect (relevant for relative astrometry as in our 
case) over a field of view of 20" is always well below 10 /xas 
and would even mainly be absorbed in the linear terms of 
the transformations. 

• Differential aberration. The classical light aberra- 
tion yields for a small field of view effectively a change in 
plate scale in one direction ( Lindegren fc Bastian|200^ I . The 
effect will be absorbed into the parameters of the transfor- 
mation when converting pixel positions to astrometric po- 
sitions via reference stars if and only if one allows also for 
the linear shear terms in the transformation. Otherwise, a 
positional error of 



6poi - cos '5 
c 



(13) 



can occur, where v ~ 30km/s is the speed of Earth and 'I' 
the angle between the apex point and the target. Numer- 
ically, this evaluates to 1 mas x cos for a field of interest 
of ±2" and hence would be a huge effect, ft is, however, 
absorbed into the shear terms of the transformation. 
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Figure 14. Summary of error sources as a function of stellar 
magnitude. The blue, horizontal lines are constant contributions. 
From top to bottom these are: residual image distortions, PSF un- 
certainty, coordinate transformations, differential chromatic aber- 
ration, differentail tilt jitter and extinction uncertainty. The red, 
diagonal, solid line is the SNR limit. The dot-dashed curves are 
for the central cusp (r = 0.5") and show from bottom to top the 
confusion error, the halo error and in thick the sum of the errors 
for that radius. The dashed curved are the same for r = 3.5". 



5 SUMMARY & CONCLUSION 

We have analyzed a multitude of error sources that poten- 
tially infiuence and bias stellar positions as obtained from 
adaptive optics assisted imaging data in crowded stellar 
fields. The domain of application in this work is the stellar 
cusp in the Galactic Center. Figure [l4] summarizes the vari- 
ous error sources for a mosaic of 48 frames with D1T= 17.2 s 
and NDIT=2 using NACO in K-band and the 13mas/pix 
image scale, assuming the current instrumental setup and 
good observing conditions. 

We find that for stars fainter than « 15 the main er- 
ror source for position measurements is halo noise, resulting 
from the imperfectly subtracted PSF seeing halos of sur- 
rounding, brighter stars. For stars of rriK ~ 14 and brighter, 
the main limitation is residual image distortions and the un- 
certainty in the PSF. In order of importance for a star with 
rriK = 15 in the central arcsecond, we estimate the following 
error contributions: 

• Halo noise: 300 /.tas 

• Residual image distortions: 130 /las 

• Confusion noise: 100 /las 

• PSF uncertainty: 50 ^as 

• SNR induced position uncertainty: 45 /ias 

• Coordinate transformations: 30 /xas 

• Differentail tilt jitter: 15/xas 

• Extinction uncertainty: 10 /^as 

• Differential chromatic aberration: 10/xas 

• Detector non-linearity: < 5 /las 

It is worth noting that our figure [14] cannot be com- 



pared to figure 3 of Chez et al. (20081, who only show the 



centroiding uncertainties and not the astrometric accuracy. 
The proper comparison is with our figure [4] 

Given the dominance of the halo noise in the error bud- 
get, the most promising route to improve the astrometry on 
existing data is to get better PSF estimates, in particular 
for the large seeing halo. This task is complicated due to 



the stellar crowding which is so dense that essentially no 
star has an unperturbed halo. It may also be worthwhile 
to think about blind deconvolution, which extracts the PSF 
while finding stars in an iterative way and in principle could 
use as much information as possible from any image. The 
advantage would be that basically all stars are used to es- 
timate the PSF halo, and not just the ones which a user 
has selected a priori. Another way to improve the PSF esti- 
mates might be the techniques of PSF reconstruction from 
the wavefront sensor data ( Gendron et al.|[2006 l. 

For the brightest stars (jtik ~ 14), the next step of im- 
provement would come from a better handling of residual 
image distortions. One possible route could be to actually 
map the distortions and then apply a correction to the pixel 
positions obtained from the data. For that, it is also impor- 
tant to have a well designed, well aligned and stable optical 
system. 

Other, desirable improvements which could be imple- 
mented in the analysis using current data are: a) take into 
account anisoplanatism during the deconvolution process 
(moderately important), b) correct for the atmospheric chro- 
matic bias due to the different stellar types (important in 
H-band only) and c) extend the linearity correction to earlier 
data obtained with another detector (barely important). 

It should also be noted that most of the errors effectively 
act as random errors, if many epochs are considered. Hence, 
we expect that observing at many epochs will average out 
many errors discussed. Nevertheless some care needs to be 
taken that this actually happens. The halo noise might get 
correlated between epochs if consistently the same PSF stars 
and same PSF size get used. Using different bands, experi- 
encing variable atmospheric conditions, the motion of stars 
and varying the way to estimate the PSF during the analysis 
help to avoid that the error turns into a bias. For the analy- 
sis, this means that from a given data set frames should be 
added up until the faintest stars of interest are well detected, 
but not beyond that. If more frames are available, it is bet- 
ter to create a second mosaic, thus lowering the halo noise 
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induced error. The only source of error which unavoidably 
leads to correlated position errors is unrecognized confusion 
events, since for typical proper motions, stars are confused 
for a few years. 

Our results also indicate that further instrumental im- 
provements mainly would come from better angular resolu- 
tion and higher Strehl ratios. It is somewhat trivially clear 
that this would yield better astrometry. The value of our 
analysis is mainly to show that this indeed is limiting the 
astrometry in the Galactic Center. Current observations are 
not SNR limited, for example. Also, we cannot exclude that 
at the level of ~ 200 /xas the astrometry has reached a limit, 
given current telescopes and adaptive optics systems. 

Since the resolution is essentially set by the apertures 
of the telescopes, only extremely large telescopes (with 30 m 
to 40 m aperture) or interferometric means in the future will 
further improve the angular resolution. 

The Strehl ratio deserves a bit more discussion. Clearly, 
the astrometry in the Galactic Center would benefit from an 
adaptive optics system that yields the highest Strehl ratio 
possible over a relatively small field. This means that adap- 
tive optics systems that yield a more moderate correction 
over wider fields are not well-suited for studying stellar dy- 
namics close to the massive black hole in our Milky Way. In- 
teresting for the GC seem the concepts for extreme adaptive 
optics systems that should deliver Strehl ratios up to 90% 
over a limited region ([Macintosh et al~ 2003 Conan et al 



20041 Finally, good astrometry will continue to depend on 
suitable atmospheric conditions. 

In a nutshell: Further improving the astonishing mea- 
surements of stellar orbits in the GC requires larger tele- 
scopes, higher Strehl ratios, better knowledge or reconstruc- 
tion of the seeing halo of the PSF, correction of residual 
image distortions - and ideally even a combination of these. 
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